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Abstract: Based on the standard Gaisser’s formula, a modified parametrization for the sea-level cosmic-ray muon 
flux is introduced. The modification is verified against experimental results. The average vertical cosmic-ray muon 
intensity as a function of depth of standard rock is simulated using the modified formula as input to the MUSIC 
code. The calculated muon intensities is consistent with the experimental measurements. 
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1 Introduction 

In 1990, Gaisser introduced a formula for describing 
the cosmic-ray muon flux at sea-level 0 i, 

0.054 

1.1 COS0 

850 GeV 

a) 

where is the differential flux in units of cm ^sr 
is the muon energy in GeV, and 9 is the zenith angle. 
This formula can be used for calculating muon-induced 
background for various underground experiments, such 
as neutrino experiments, double-,5 decay experiments, 
dark matter search experiments, and the like. However, 
this standard Gaisser’s formula is only valid under the 
following two conditions 0: 

1. the curvature of the Earth can be neglected {6 < 
70°). 

2. muon decay is negligible > 100/cos9 GeV). 

In this note, we suggest a modified parametrization 
that can overcome these shortcomings. 

2 Modified parametrization 

As stated in condition 1 in the introduction, the stan¬ 
dard Gaisser’s formula cannot be applied to all zenith 


angles. According to references 0[1|! when taking the 
curvature of the Earth into account, the observed zenith 
angle on the ground, 9, and the zenith angle at the pro¬ 
duction point of the muon in the atmosphere, 9*, are 
different. Figured] illustrates the relation between these 
two angles. 



Fig. 1. The relation of the observed zenith angle 
of muons, 9 *, to the zenith angle at the muon 
production point in the atmosphere, 6. R is the 
radius of the Earth. Adopted from [3[J] 

Chirkin (originally tabulated in Volkova (1969)) 0 
pointed out that the quantity cos(0*) can be related to 
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COS0 such that 


cosO* = 


(cos0)^ + + P2 (cos0)^3 +P4(cos0)^5 


1+P?+P2+P4 


( 2 ) 

where Pi, P 2 , P 3 , Pa, and P 5 are parameters given in Ta- 
ble[TJ Equation[2]is indeed a convenient parametrization 
of the effect of the Earth curvature. 


Table 1. Parameters in Eq. laa. 


Pi P2 

Pa 

P4 

Ps 

0.102573 -0.068287 

0.958633 

0.0407253 

0.817285 


As shown in Figure [2l due to the seconds reason 
stated above, the standard Gaisser’s formula cannot de¬ 
scribe the experimental results well at low energies. We 
have modified Equation [T] by adding a term to the 
factor. When the muon energy increases, this term 
becomes negligible and the original functional form of 
Gaisser is recovered. To determine the parameters of the 
new term, we fit the modified equation with the zenith 
angle given in ^nation [3 to the world cosmic-ray muon 
measurements i0 i 1 0. The resulting modifica¬ 
tion is Equation 
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level which then served as the input to the MUSIG code 


ll| for transporting the simulated muons to a specific 


depth of standard rock. The predicted muon flux is com¬ 
pared to the experimental data of vertical muon intensity 
at different depths of rock overburden in Figure[31 where 
the data are taken from Reference 12|. The simulated 
results and experimental data agree well. 



Fig. 3. Average vertical muon Intensity versus 
depth of standard rock. Black points are exper¬ 
imental data from reference [l^. Red solid tri¬ 
angles stand for the simulated results using the 
modihed parameterization (Equation [3|. Green 
hollow triangles are the simulated results using 
the standard Gaisser’s formula (Equation [T]) . 


3 Conclusion 


where the parameters 3.64 and 1.29 are obtained from 
the fit. The description of Equation |3] is compared to 
the measured cosmic-ray muon fluxes in Figure [5J The 
modified parametrization matches the experimental re¬ 
sults fairly well for different zenith angles. 



Fig. 2. Comparison of the modified parametriza¬ 
tion to the measured cosmic-ray muon fluxes.The 
experimental data are from i 0 i i 0 


To validate the modified parametrization, Equation[3] 
was used to generate the cosmic-ray muon flux at sea 


We have obtained a modified Gaisser’s formula that 
extends the range of applicability to all zenith angles and 
lower energies. The new parametrization can be used 
conveniently and reliably for representing the cosmic-ray 
muon distribution at sea level for ground detectors, as 
well as for underground experiments after it is coupled 
to a software package for transporting the surface muons. 
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